Start populations in the lower third of the simulation and let them move as they co-evolve. It is possible that as newts and snakes move into new territory they might create more spatial phenotype correlation.
What does it mean to be correlated? Might this occur if I simulate some controlled movement?
I am running a coevolution experiment where newts and snakes interaction with each other over long periods of time. When they interact the newt or snake with the highest level of toxicity or resistance will survive. As newts and snake interact across space, they develop exaggerated traits that can become spatial correlated. My experiments examine when traits became exaggerated or correlated.
I created a simulation study to observe the co-evolutionary outcome of the newt-snake interaction with different genetic architectures (GAs) in a spatial setting. I hypothesized that we would see an interaction (co-evolutionary arms race) between newt and snake phenotype under some GA combinations when newts and snakes were evolving over geographical space. Each GA is paired with another GA creating 16 combinations.
GA1 experiment values:
Each GA combination and trial has its own msprime simulation, but the msprime file is shared between the 5 section and 7 section slim runs.
## All cor, lit, and grid files exist!
## This program will now end!
## Group.1 x
## 1 1e-08_0.005_1e-08_0.005 -0.03323921
## 2 1e-08_0.005_1e-09_0.05 -1.30511227
## 3 1e-08_0.005_1e-10_0.5 -2.50095446
## 4 1e-08_0.005_1e-11_5 -1.19918146
## 5 1e-09_0.05_1e-08_0.005 1.33051637
## 6 1e-09_0.05_1e-09_0.05 -0.43179017
## 7 1e-09_0.05_1e-10_0.5 -1.13158335
## 8 1e-09_0.05_1e-11_5 -0.48635858
## 9 1e-10_0.5_1e-08_0.005 2.03584772
## 10 1e-10_0.5_1e-09_0.05 0.39730616
## 11 1e-10_0.5_1e-10_0.5 -0.73000969
## 12 1e-10_0.5_1e-11_5 -0.40944406
## 13 1e-11_5_1e-08_0.005 0.98550774
## 14 1e-11_5_1e-09_0.05 -0.79978450
## 15 1e-11_5_1e-10_0.5 -1.42371425
## 16 1e-11_5_1e-11_5 -0.60451822
## [1] "pattern 1e-08_0.005_1e-08_0.005_3"
## [1] "Cor between average snake pheno and local cor -0.423984702204067"
## [1] "Cor between average newt pheno and local cor -0.262198745181791"
## [1] "Cor between average dif pheno and local cor 0.200962706000563"
## [1] "Cor between newt pheno and snake 0.688195150229443"
## [1] "pattern 1e-08_0.005_1e-08_0.005_1"
## [1] "Cor between average snake pheno and local cor -0.0780136619989221"
## [1] "Cor between average newt pheno and local cor -0.0304579460944814"
## [1] "Cor between average dif pheno and local cor -0.104524490191293"
## [1] "Cor between newt pheno and snake 0.885641401406511"
## [1] "pattern 1e-09_0.05_1e-09_0.05_2"
## [1] "Cor between average snake pheno and local cor 0.423145833547384"
## [1] "Cor between average newt pheno and local cor 0.594390934146489"
## [1] "Cor between average dif pheno and local cor -0.614239705996602"
## [1] "Cor between newt pheno and snake 0.85059971017444"
## [1] 0.4002541
## [1] 0.4181679
## [1] 0.1298847
## [1] 0.5833832
## [1] 0.200901
## [1] 0.1563885
## [1] 0.09847233
## [1] 0.409964